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Abstract 

The many-body space fractional quantum system is studied using the density 
matrix method. We give the new results of the Thomas-Fermi model, and obtain 
the quantum pressure of the free electron gas. We also show the validity of the 
Hohenberg-Kohn theory in the space fractional quantum mechanics and generalize 
the density functional theory to the fractional quantum mechanics. 

1 Introduction 

The famous Feynman path integral l[T]|2l reveals the relationship between frac- 
tals f3l and quantum mechanics. The background of the Feynman approach to quan- 
tum mechanics is a path integral over Brownian paths. The Brownian paths are non- 
differentiable, self-similar curves that have a fractal dimension that is different from its 
topological dimension. Recently, Laskin |4-8 1 used Levy paths instead of the Brow- 
nian ones in the path integral and obtained a space fractional Schrodinger equation 
(FSE) and then developed fractional quantum mechanics. Following him, some other 
generalizations to the standard Schrodinger equation appeared. Naber |9| showed a 
time fractional Schrodinger equation , and then Wang and Xu 1 10| combined the two 
kinds of fractional Schrodinger equations together to construct a space-time fractional 
Schrodinger equation. Moreover, the Fractional Heisenberg equation IfTTHTSl is also 
developed. These generalizations may describe more extensive fractal phenomena in 
the microcosmic world. 

There is a physical reason for the merger of fractional calculus P14','151 with quan- 
tum mechanics. The Feynman path integral formulation of quantum mechanics is based 
on a path integral over Brownian paths. In diffusion theory, this can also be done to gen- 
erate the standard diffusion equation; however, there are examples of many phenomena 
that are only properly described when non-Brownian paths are considered. When this 
is done, the resulting diffusion equation has factional derivatives llT6l - IT8l . Due to the 
strong similarity between the Schrodinger equation and the standard diffusion equation 
one might expect modifications to the Schrodinger equation generated by consider- 
ing non-Brownian paths in the path integral derivation. This gives the time-fractional, 
space-fractional, and space-time-fractional Schrodinger equation 191 011 . 
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The space-fractional Schrodinger equation (4] obtained by Laskin reads (in three 
dimensions) 

dMr, t) 

ih^^^^H^^irj), (1.1) 
at 

where if/ir, f) is the time-dependent wave function , and Ha {I < or < 2) is the fractional 
Hamiltonian operator given by 

Ha^-Daimr + V{r,t). (1.2) 

Here Da with physical dimension [Da] - [Energy]'"" x [Length]'' x [Time]"" is de- 
pendent on a [Da - 1 12m for a - 2, m denotes the mass of a particle] and (^V)" is the 
quantum Riesz fractional operator ISJUS) defined by 

(^V)"(/r(r, f) = ^—T f dVe'P'^^lpr r e-'P'^VCr, Od^r. (1.3) 

(2niiy J J 

So far, the research on the fractional quantum mechanics is based on solving the 
FSE. In Ref. ||4l[T9H2T1. the exact solutions to the FSE with free particle, Dirac delta 
potential, linear potential, and some other simple potential fields are given. However, 
because of the double complexity of the fractional calculus and quantum mechanics, 
it is hard to solve the FSE directly for most cases, especially for the many-body quan- 
tum systems. In this paper, we use the density matrix method to study the many-body 
fractional quantum system. We give the density matrix description of the fractional 
quantum system, study the Thomas-Fermi model under the framework of the fractional 
quantum mechanics, and generalize the density functional theory to the fractional quan- 
tum mechanics. 



2 Density matrix description of the fractional quantum 
system 

We are interested in systems of identical particles (electrons, say) moving in a 
given external field and interacting with each other with pair forces. In most cases, 
one is concerned with atoms and molecules without time-dependent interactions, so 
we may focus on the time-independent FSE |T9l |. For an isolated N-electron atomic 
or molecular system in the Born-Oppenheimernon-relativistic approximation Il22ll23l . 
the FSE is given by 

HaW = E\\i, (2.4) 

where E is the electronic energy, \\i = y(xi,X2, ■ ■ ■ ,xpf) is the time-dependent wave 
function ( the coordinates x, of electron / comprise space coordinates r, and spin coor- 
dinates s, ), and Ha is the fractional Hamiltonian operator, 

Ha^t + V„e + V,e, (2.5) 

where t = Dai-h^Ai)"^^ is the kinetic energy operator, V„e = 2/1 1 '-'(r/Ois the 

electron-nucleus attraction energy operator, in which ^(r,) = - 2i> is the external 

rm 

potential acting on electron /, the potential due to nuclei of charges Z,j, and Vee = 
— is the electron-electron repulsion energy operator. The total energy W of this 
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system is the electronic energy E plus the nucleus-nucleus repulsion energy, V„„ - 
T.»<e That is, W = £ + y„„. 

Now we introduce the density matrix Il22ll23l to describe this quantum system. 
Considering the electrons system in a pure state YAr(xiX2 ■ ■ ■ x^r), then we can define 
the density matrix as 

yA'Cx'iXj ■ ■ • x^, X1X2 ■ • ■ Xiv) = Yiv(xix2 ■ ■ • xJy)\|;^(xiX2 • ■ ■ x^v), (2.6) 

which can be viewed as the coordinate representation of the density operator, 

JN = IViyXVA^I, (2.7) 

where the Dirac notation is used. Note that yf^ is a projection operator, and for normal- 
ized M/m, 

ti-(fw) = J W(x'^)v;^(x^)flfx^ = 1, (2.8) 

in which tr(7A'), the trace of the operator jm, is defined as the sum of diagonal elements 
of the matrix representing yw, or the integral if the representation is continuous. In the 
state the expectation values of any observable A can be calculated by ll22l 

(A) = tr(ArM) = tr(rivA), (2.9) 

where A is the Hermitian linear operator for the observable A. In view of Eq. ( |2.9l l, the 
density operator jj^i carries the same information as the N-electron wave function |\|/). 
yjv is an operator in the same space as the vector Note that while is defined only 
up to an arbitrary phase factor, for a state is unique, also is Hermitian. Further 
more, we can define the reduced density matrices for fermion systems ll22l . 



7;,(x'iX^ • ■ ■ x;,, X1X2 ■ ■ ■ ''^^ = (p) J" ■ ■ ■ J" rA'(x'iX^ 

Xp+1 ■ ■ ■ Xa,, X1X2 ■ ■ ■ X,, ■ ■ ■ XN)dXp+i ■ ■ ■ dXN, 



(2.10) 



where is a binomial coefficient. In particular. 



72(x'iX2,XiX2) = ^^^^ J ■■ ■ J yNix[X2Xi ■ ■ ■Xn,XiX2X2 ■ ■ ■ Xiv)<iX3 ■ ■ ■ dx^. 



(2.11) 



and 



yi(x'i,xi) = J^--- J"yiv(XjX2---Xiv,XiX2---Xiv)c/x2---<iXiv. (2.12) 

Many operators of interest do not involve spin coordinates, for instance the Hamil- 
tonian operators for atoms or molecules. This makes desirable further reduction of 
the density matrices of the first and the second orders. We define the first-order and 
second-order spinless density matrices by 

Pi(r'i,ri) = J"yi(r'iSi,riSi)t/si. (2.13) 

and 

P2(r'ir2,rir2) = ^r2(r'iSir2S2,riSir2S2)c/sic/s2. (2.14) 
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We also introduce a shorthand notation for the diagonal element of pi and pi'- 

p(ri) =pi(ri,ri),p2(ri,r2) = P2(rir2,rir2). (2.15) 
Note that the diagonal element of p\ is just the electron density, because that 

p(ri) =pi(ri,ri) = ^ ■ ■ ■ ^ \y]f\^dsidiL2- ■ ■ diiN- (2.16) 

A useful relationship between the first and the second order spinless density matrices 
is given as follows, 

Pi(r'i,ri) = -^^—^ J p2(r',r2,rir2)t/r2. (2.17) 

In particular, 

P(ri) = J P2(ri,r2)t/r2. (2.18) 

After these definitions, the expectation value, for an antisymmetric N-body wave func- 
tion, of a local one-electron operator Oi - YJiL\ Oi(ri), can be expressed as 

<Oi) = tr(()if^) = J[Oi(ri)p,(r;,ri)]r;=r,^ri. (2.19) 

Similarly, the expectation value of a two-electron operator 02 = Yj^kj Oiiri, rj) is 

{62} ^ HOiYn) ^ j [02(rir2)p2(r'ir2,r,r2)]r;=r,,r-r2t^rit/r2. (2.20) 

Using the above results, combining all the parts, the expectation value of the Hamilto- 
nian i2.5\ . is obtained, 

E =tr(H„7A,) = £'(pi(r'i,ri),p2(ri,r2)) 

= J[D,(-h^Arr'Wr',r)]r,=rdr + J v(r)p(r)dr+ JJ ^^p2{ri,r2)dr,dr2. 

(2.21) 

The three terms in this formula represent respectively the electronic kinetic energy, 
the nuclear-electron potential energy, and the electron-electron potential energy. If we 
write 

P2(ri,r2) = ip(ri)p(r2)[l + h{rur2)]. (2.22) 
Then the third term in Eq. ( 12.211 1 can be rewritten as 

~P2(ri, r2)iir2iiri = 7[p] + ^ ^ ^p(ri)p.„.(ri, r2)iiri Jr2, (2.23) 

where ^ 

J{p]-\^^ ^^P(ri)p(r2)dridr2, (2.24) 

and Pxc(ri,r2) = p(r2)/i(ri,r2). Here, Pxc(T^i,r2) is called the exchange-correlation 
hole 11^123], and 



P,,(ri,r2)t/r2 = -1. (2.25) 



So far, the density matrix description of the fractional quantum mechanics has been 
made. Now we begin to study the fractional quantum mechanics by use of these results. 
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3 The Thomas-Fermi model 



In the standard quantum mechanics, the first attempt to study the density functional 
theory of electron structure is the works of Thomas and Fermi in the 1920s. Here, under 
the framework of the fractional quantum mechanics, we also start from the Thomas- 
Fermi Model Il22l424l . We divide the space into many small cells, each of side I and 
volume Ay = P, each containing some fixed number of electrons AA^, and we assume 
that the electron in each cell behave like independent fermions at the temperature OK, 
with the cells independent of one another. Then we first study the energy levels of a 
particle in each cell. The particle is like in a three-dimensional infinite well, 

f 0, ifx,3;,ze (0,0; 
V(x,y,z)^\ ^ . (3.26) 

oo, otherwise. 

So we have the following FSE, 

Da(-h^Ar'\(,r) = £v(r), (3.27) 

with the boundary conditions 

\|/(r) = 0, when x = /, or y = /, or z = /, or x = y = z = 0. (3.28) 

It is easy to prove that \\i(r) = Ce''' '" is a basic solution to Eq. ( 13.271 ). Here, k = 
(kx; ky, k.), satisfying 

E = D^h^M'. (3.29) 
The general solution to Eq. (13.271 1 can be expressed by 

\|/(r) = A sin(k ■ r) + B cos(k ■ r). (3.30) 

Considering the boundary condition ( 13.281 1, we have 

B = 0, sin(;t^L) = sin(;t,.L) = sin(;t,-L) = 0, (3.31) 

which means 

kxL - rijcTT, kyL - riyTt, k^L - n,n, (3.32) 

n 

where Hi, n,,, Mj. can be any integer So k = {kx,ky,kr) - —(nx,ny,n.). Then from 
Eq. (Il29] l. we get 

£(«„ riy, «,) = D,X\kr = D,{—rR\ (3.33) 

where R - ^n]. + + n^. Considering the spin states of electrons, every space state 
(n V, ny, Hj) can be occupied by two electrons. For high quantum numbers, that is, for 
large R, the number of distinct energy levels with energy smaller than E{n y, «v, n,) can 
be approximated by the volume of one octant of a sphere with radius R in the space 
(n I, Hy, n~). That's because that every array of positive integers {rix, riy, corresponds 
to one point in the first quadrant of the space (n^, n-) and there is only one grip point 
in per volume unit. Therefore, the number of distinct quantum states, in the interval 
{R, R + dR), is given by 

dN = 2-]-- AnR^dR = nR^dR, (3.34) 
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where the factor 2 enters because each energy level is doubly occupied,by one electron 
with spin a and another with spin /3. By the Pauli exclusion principle, we know that 
the number of distinct quantum states is also the number of the electrons. When is 
large, the change of can be viewed as continuous. Differentiating Eq. (13.331 1 yields 

dE = aDai^)"R"-^dR. (3.35) 

Thus 

dN ^nR^aD^y^—TR^-^dE, (3.36) 
nn 

That is 

^^nR\aD,,)-U-^rR'-\ (3.37) 
dE nn 

the right side of which gives the number of electrons with energy in the interval (E, E + 

dE). Thus, the number of electrons with energy smaller than the Fermi energy Ef is 

= f ' ^dE = l(D,rV"(-Lf(Ef)y". (3.38) 
Jo dE 3 nn 

The total energy for the volume cell can now be calculated as 

^ = r E^dE = ^(Dar'^"{^)\Efr^''^l\ (3.39) 
Jo dE a + i nn 

which can be rewritten as 

E = -^NEf = -l-D,n"{3n^rl\^r^^^'L\ (3.40) 
ff + 3 a + 3 

N N 

When the volume cell V — > 0, = dV - dr and electronic density n - — = — = 

, , ^ ^ V 

p(r). From Eq. ( 13.401 ). we obtain the total energy dE for the volume cell dr as 

dE = -^DaW{2>n^fl^(p{Y)fl^^''dr. (3.41) 
a + 3 

Adding the contributions from all cells, we find the total kinetic energy to be 

TtfIp] ^CfJ (p{r)r''^'dr, (3.42) 

3 

where Cf DaK'i'in^)"^^ . This is the kinetic energy functional for the Thomas- 

a + 3 

Fermi model in the fractional quantum mechanics. When a = 2, by use of the atom 
unit, we have Ttf[p\ = Cf J(fi(r))^^^dr, and Cf = ^(371^)2/3 ^ 2.871, which accords 
with the standard quantum mechanics Il22l . 

Recalling Eq. ( 12.21b . we can get an energy formula for an atom in terms of electron 
density alone, 

ETF[p(r)]^CF J(p{r)r''^'dr-Ze^ J + jj ^^p2{rur2)drydr2. (3.43) 

For the ground state of an atom of interest the electron density minimizes the energy 
functional £'j-/r[p(r)], under the constraint 



N = N[p{r)] = p{r)dr, (3.44) 
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where N is the total number of electrons in the atom. By the method of Lagrange 
multipliers, the ground-state electron density must satisfy the variational principle 



6 ^EjFlpir)] - fiTF (J P(r)flfr - Njj = 0, (3.45) 
where htf is the Lagrange multiplier. After executing the variation, we can get 

f^TF = — Cfp'^^^(r) - ^(r), (3.46) 



where (p{r) is the electrostatic potential at point r due to the nucleus and the entire 
electron distribution: 

^(r) = r^^'"^- (3.47) 

r J |r-r2| 

Eq. (13.461 1 can be solved in conjunction with the constraint (13.441 ). and the resulting 
electron density then inserted in ( 13.43b to give the total energy. There are the results 
of the Thomas-Fermi model for the atom in the fractional quantum mechanics. Just as 
in the standard quantum mechanics, it is a simple model, without considering the ex- 
change and correlation energy terms. But in all the results, the parameter a is contained. 
Proper choose of a may provide better results than the standard quantum mechanics. 

Now we turn to solve the equations ( 13.461 1 and ( 13.471 ). For neutral atoms, when 
r — > oo, i^(r) — > and p(r) — > 0, so from Eq. ( 13.461 ). we have fj.TF = 0, which means 

p(r) = C-^'"if^'"(r), (3.48) 

where C = Cf (a -i- 3)/a. Then using the results of the Poisson equation Il25l . the 
expression of (f{r) in Eq. ( 13.47b satisfies the following equation 

VV(r) = Ane^pir) - Ane^Zdir), (3.49) 

Substituting Eq. ( 13.48b for p(r) gives 

VV(r) = 4-ne^C-^^''ip^'''(r) - Ane^ZSir). (3.50) 

Additionally, by the spherically symmetric of the atom, we can write 

ip{r)^^{r)^—oj{r). (3.51) 
r 

Then we have 

o Ze^ d^cj(r) 
VV(r) = - Ane^Z(jj{0)5{r). (3.52) 

Eq. (13.50b minus Eq. ( 13.521 ) side by side gives 

—11 ^ Mr'-nc^{r)]-^ , (3.53) 

where 

M ^AnC-^I^Z^^^-^e^l", 

and the boundary condition w(0) = 1 is introduced. Letting x - ar and a - M"^'^*" 
we obtain a differential equation for w(x), 

d^(jL>(x) , 3 3 

^^=xi-«Mx)]«, (3.54) 
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Figure 1: The curves for the numerical solutions of a;(x) for a - 1.1, 1.5,2. a - 2 
corresponds to the standard quantum mechanics (see Ref. ll22^B4l '). 



with the boundary conditions oj(0) — 1 , w(oo) - 0. This problem can be solved 
numerically. The numerical result is given in Fig. [T| for different a. From the figure, 
we know that for all a, the curves decrease monotonously, and for a smaller a, the 
curve decreases more rapidly at the beginning, but more slowly as x — » oo. Note 
that Eq. (13.54b has an exact solution (l)(x) - Bx^, where A - 3(a - l)/ia - 3), and 
B - [A{A - This solution does not satisfy the boundary condition at x = 0, 

but can be a good approximate expression for oj{x) when x is large. 

4 Quantum pressure of the free electron gas 

By use of the results in the previous section, we can study the quantum pressure 
of the electron gas in the fractional quantum mechanics. Defining the variation of the 
volume of the electron gas by dQ., the needed energy from outside by dA, then the 
pressure p, of the electron gas, satisfies 



dD. = -pdA, 



(4.55) 



Meantime, the internal thermal energy increases by dU — dA, so 



dU 



(4.56) 




For degenerate Fermi electron gas, using Eq. ( 13.391 ), we have 



3N 



(4.57) 



U = NE = 



a + 3 



Thus 



dU = 



3N 



■dEf. 



(4.58) 



a + 3 



Remembering Eq. ( 13.381 1. we can get 



1 3 3 

E" = N-iD,)-' 

■I TT 



1 {nnf 



(4.59) 



Q 
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where the volume of the electron gas O = L . The above equation can be converted 
into 

3 

— \nEf - constant -\nQ,. (4.60) 

a 

Differentiating the both sides gives 

3 Q. 

dQ. = dEf. (4.61) 

aEf 

Then, with the help of Eqs. (14.58b and (14.611 1, the pressure p of the electron gas can be 
calculated by Eq. ( I4.56l i, 

dU a N a 



Ef= pEf, (4.62) 

where p = N/D. denotes the density of the electron gas. Finally, combing Eq. (14.59b 
and ( 14.62b , the pressure p is expressed by the density p as follows, 

(3;r^)?Da.rpt+'. (4.63) 



a + 3 



This means there is a stabilizing internal pressure in the solid object. But this pressure 
derives ultimately from the anti-symmetrization requirement for the wave functions of 
identical fermions, so it can be called the degeneracy pressure |26|. When a = 2, 
Eq. (14.63b reduces to the result in the standard quantum mechanics 1.26 J . 



5 Accurate theory: Hobenberg-Kohn theorems 

In the previous sections, we study the many-body fractional quantum system under 
the Thomas-Fermi model, an exquisitely simple model. To get the accurate results, 
we recognize the density functional theory Il22ll23]| just as in the standard quantum 
mechanics. From the discussion in Ref. (4\, we know that the hamiltonian for the 
space fractional quantum system is still a hermite operator Employing the minimum- 
energy principle for the ground state, we can easily prove that the Hohenberg-Kohn 
theorems |22,23 | still hold in the space fractional quantum mechanics. We now give 
the two theorems without proof, and the proof is exactly the same with the standard 
quantum mechanics: 

1 .The electron density p(r) determines the ground-state wave function and all the other 
electronic properties of the system. 

2. For a trial density p(r), such that p(r) > and J p(r)dr - N, 

Eo < E,[p], (5.64) 

where Ey[p] is the energy functional of the quantum system with external potential 
u(r). From Eq. ( I2.21l )in Sec.|2] the total energy can be written as 

E,[p]^T[p] + V„e[p] + Vee{p]^ j v(r)p(r)dr + F hkIpI (5.65) 

where FhkIp] - T[p] + yee[p].We may write y(.f.[p] = J[p] + nonclassical term, where 
J[p] is the classical repulsion. The nonclassical term is a very elusive, very impor- 
tant quantity; it is the major part of the "exchange-correlation energy". Assuming 
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differentiability of Ey[p], the ground-state electron density must satisfy the variational 
principle 



6 |£Jp(r)] -fi{^j p(r)dv - ^ - 0, (5.66) 

which gives the Euler-Lagrange equation 

SEt, 6Fhk\p\ 
II = ^-r-: = v(r) + ^ . . ■ (5.67) 
(5p(r) 5p{r) 

The quantity yU is the Lagrange multiplier and called the chemical potential. If we knew 
the exact Fhk[p\, Eq- ( 15.66b would be an exact equation for the ground-state electron 
density. As in the standard quantum mechanics, accurate calculational implementations 
of the density functional theory are far from easy to achieve, because of the unfortunate 
fact that the functional Fhk[p\ is hard to come by in explicit form. To implement 
density functional theory, the Kohn-sham method 1122112311 can also be used here and 
the resulting Kohn-sham equation will contain a fractional operator. Further research 
on this subject will be done in our future work. 



6 Conclusions 

In this paper, the many-body space fractional quantum system, governing by frac- 
tional Schrodinger equation obtained by Laskin, is studied using the density matrix 
method. We give the results of the Thomas-Fermi model, calculate the quantum pres- 
sure of electron gas in the fractional quantum mechanics. We also show the validity 
of the Hohenberg-Kohn theorems, and generalize the density functional theory (DFT) 
to fractional quantum mechanics. These results provide the possibility to study the 
fractional quantum system without solving the FSE directly. 

The density functional theory is a popular tool applied to study the electron struc- 
ture of solids Il27l - |29l . but because the exchange-correlation term is unknown and hard 
to define, there are always modifications needed. The fractional quantum mechanics, 
which is built in the sense of Levy flight, is a generlisation to the standard one, which 
is built in the sense of Brownian motion. The Levy flight is a natural generalization 
of the Brownian motion, and the fractal dimension of Brownian paths is two, while 
that of Levy paths is a (1 < a < 2). The many-body quantum systems can exhibit 
the characteristic of fractal with the fractal dimension not equal to two, so it is reason- 
able to believe that the fractional Schrodinger equation can better describe the behavior 
of the microcosmic particles. The density functional theory in the fractional quantum 
mechanics may be an alternative way to study the structure of the many-body quantum 
systems, and without making more modifications to the exchange-correlation terms, by 
choosing a proper parameter a, the DFT in the fractional quantum mechanics may pro- 
vides better results than in the standard quantum mechanics, which needs our further 
studies. 
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